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Abstract 

The question about the existence of Bogoliubov's quasiparticles in the BCS wave functions underneath Gutzwiller's projection is 
of importance to strongly correlated systems. We develop a method to examine the two-particle excitations of Gutz wilier-projected 
BCS wave functions by using the variational Monte Carlo approach. We find that the exact Gutzwiller-projected quasiparticle 
(GQP) dispersions are quantitatively reproduced by the Gutzwiller-projected Bogoliubov quasiparticles (GBQP) except the regions 
where cZ-wave Cooper pairing is strong. We believe GBQP provides a reasonable description to the low-energy excitations in 
strongly correlated superconducting systems because GBQP becomes more stable than GQP near the antinodes. In addition, the 
intimate connection between Gutzwiller's projection and af-wave Cooper pairing may also imply that strong correlations play a 
significant role in the nodal-antinodal dichotomy seen by photoemission experiments in cuprates. 

Keywords: A. superconductors; C. variational Monte Carlo; D. strong correlation; D. quasiparticle excitations 
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1. Introduction 

Two of the most intriguing puzzles in the study of high T c 
cuprates are the unexpected non-BCS behavior and the non- 
quasiparticle nature in the superconducting states and the nor- 
mal states, respectively 01 • To resolve those puzzles, the rele- 
vant low energy physics based on the projection out of the de- 
grees of freedom at high energy must be embedded in a doped 
Mott insulator [2, 3]. The Gutzwiller-projected BCS wave func- 
tion is the appropriate description of the superconducting state 
in cuprates |4|], while strong correlations make theoretical ap- 
proaches extremely difficult. However, based on the frame- 
work of the Gutzwiller-projected states, the issues related to 
the finite-temperature physics of cuprates are still unclear |5j|6(]. 
Therefore, the first step in studying the excitations is to under- 
stand the structure of low-lying quasiparticle excitations. 

It has been experimentally observed that the low-lying ex- 
citations of superconducting cuprates resemble BCS Bogoli- 
ubov's quasiparticles (BQP) [7]. Also, many theoretical stud- 
ies on the Gutzwiller-projected BQP (GBQP) excitations have 
been presented few years ago fi U E G2 [HI Q • The 
ansatz used for the GBQP excited states is based on the suc- 
cess of the Gutzwiller-projected BCS wave function. Owing to 
exact diagonalization results indicating the well defined BCS- 
like BQP as low-energy excitations of the t - J model [8], the 
GBQP excited state given in Eq.© is expected to be a simply 
renormalized BQP excitation despite lack of analytical proof. 
Relied on the careful fitting simulations foil , we found they 
are quantitatively satisfied with the renormalized BQP picture. 
Even so, we still have no knowledge of the exact Gutzwiller- 
projected quasiparticles (GQP). On the other hand, to explain 
some unusual features seen by angle-resolved photoemission 



spectroscopy (ARPES) in cuprates, the extension beyond the 
single-mode approach has also been studied I115h . Therefore, 
from a theoretical point of view, the difference between GBQP 
excitations and GQP excitations should be clarified. 

Let us briefly summarize the key messages involved in this 
article. We begin by detailing the procedure that constructs 
the two-particle excitations by using the usual GBQP picture 
and the GQP excitations in the variational Monte Carlo (VMC) 
calculation. We demonstrate the projected two-particle excita- 
tion is reasonable to be constructed by applying BQP operators 
to strongly correlated superconducting ground states. It is no- 
ticed that there is the discrepancy between GQP and GBQP near 
antinodal regions. This discrepancy arising from a close rela- 
tion between Gutzwiller's projection and c/-wave pairing may 
provide a clue to the causes of the nodal-antinodal dichotomy 
observed by ARPES measurements Ill6lll7ll . 



2. Theory 

Let us begin by 
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where the hopping fy = t, t', and t" for sites i and j being the 
nearest, second-nearest, and third-nearest neighbors, respectively. 
Other notations are standard. We restrict the electron creation 
operators c? to the subspace without doubly-occupied sites. In 
the following, the bare parameters in the Hamiltonian are set to 
be (t',t",J)/t = (-0.3,0.15,0.3). Two holes are doped into the 
extended t — f — t" — J model in 16 x 16 lattice. 
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The well-known candidate for the ground state HE 
I20I1 is the of-wave resonating-valence-bond (d-RVB) wave func- 
tion with Jastrow factors of the form 



I<Dq> = PnPjPg 



'''k + V k4 T C -k|)l >' 



(2) 



where the coefficients and are the BCS coherence fac- 
tors. The trial wave function has three projections: P^ e to 
fix the number of electrons N e , the Gutzwiller projector Pc(- 

(1 - «;T n 4)) to enforce no-doubly-occupied sites, and charge- 
charge Jastrow factors Pj to repel neighboring holes (see the 
details in Ref. 12011 '). It is not shown that the conclusions in this 
work would not be changed by P j, and hence we will ignore 
the charge-charge Jastrow factors in the following. 

Even so, it is still not easy to construct the single-particle 
excited state of the extended t - J Hamiltonian in the canonical 
ensemble. Based on Eq.(f2]), however, the simplest way is to 
define a single-particle excitation under Gutzwiller's projection 



l*hr> = Pk±iPg7{JBCS), 



(3) 



where y kcr (- u kC kcr - crv k c- k a-) is the creation of the BQP and 
cr spin index (cr = -<f). In what follows, to avoid the confusion 
due to mixing Hilbert space of different particle numbers, we 
introduce the partial particle-hole transformation to change the 
representation from (c) to (df) 112 1 L l22fl : 



(c) 



(df) 

fl 



(4) 



Two different particles, d and /, are thus introduced instead of 
down- and up-spin electrons (see the details in Ref.[22]). We 
start from the wave function without Gutzwiller's projection. 
First, the BCS wave function can be transformed into the repre- 
sentation (df), 



\BCS) -> W (u k f! k - v k d f _ k ) \0) (df) , 



(5) 



where the subscripts indicate different representations. Then 
the single-particle BCS excited state is similarly transform into 



y^BCS) -» Y] KA - V- q ) I0>W), 

q^k 



yl H \BCS > -> dlj\ Y] (u q fl q - v q <£ q ) \Q) m . 
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(6) 



(7) 



In the representation (df), the total particle number of the single- 
particle BCS excitation is no longer confused with N e + 1 or 
N e - 1 like Eq.© in the representation (c), but fixed to N - 1 in 
Eq.© and JV+ 1 inEq.©. 

On the other hand, we adopt the similar route to write down 
the two-particle BCS excitation y k ^y_ ki \BCS) in the represen- 
tation (df): 



{u k d\ + V k /j k ) Y\ ("q/iq - V-q) l°>W> 
q#k 



(8) 



According to the above equation, if we define two states |T^) 
and as 

q^k 

= /I k nMl-Vtq)|0>(^ (9) 
q#k 

the BCS ground state and the two-particle BCS excitation will 
be obviously given by 



\BCS) 



v k |^k> + M kl^{), 



yllAl lBCS) = «kl^) - v k |T{), 



(10) 



respectively. Owing to (T^T*^) = 0, Eqs.dTOt simply represent 
that both the BCS ground state and the first BCS excited state 
are able to be expanded by the two orthonormal states, l*!*^) and 

Applying the similar idea in the projected case, we can write 
down the d-RVB ground state and the corresponding first GQP 
excited state |<I> k ) (shown in Ea.(fT5b) by using the following 
two states: |<I> k ) and |<I> k ) given by 



(11) 



As well, we can continue the single-particle excitation shown 
in Eq.© to create two-particle GBQP excited state: 
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PNPcy^y^BCS). 



(12) 



Some details in the VMC calculation should be noticed. To 
avoid the divergence from the nodes in the trial wave functions, 
the boundary condition we use is the anti-periodic boundary 
condition along both x and y directions. In order to achieve a 
reasonable acceptance ratio, the simulations consist of a com- 
bination of one-particle moves and two-particle moves. Varia- 
tional parameters in the d-RVB state are optimized by using the 
stochastic reconfiguration method [23]. All physical quantities 
are calculated with the optimized parameters. We also take a 
sufficient number of samples (= 2 x 10 5 ) to reduce the statisti- 
cal errors, and keep the sampling interval (~ 40) long enough 
to ensure statistical independence between samples. 

3. Results and discussion 

Since the projections are included in Eqs. tfTTt . the trial states 
|<t> k ) and \(£> k ) are no longer orthogonal. We need to diagonalize 
a 2x2 Hamiltonian matrix in the subspace spanned by Eqs.dTTTi. 
In principle, we can reconstruct two orthogonal states, |<I> k ) and 

A f 

|<t> k ), for each momentum by using Gram-Schmidt method. The 
Hamiltonian matrix in this subspace is given by 
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df 
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where fly = (<I> k |fl|<i> k } and i,j = d or /. We further diago- 
nalize H m b to obtain the eigenstates |O k ) and |O k ) as a linear 
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Figure 1 : (a) Two-particle dispersions of different excited states for the t — t 1 — 
t" -J model doping two holes in 16x 16 lattice. The orange dashed line denotes 
the rf-RVB state (E dRVB = -90.93/). The symbols for momenta: T = (0.5, 0.5), 
X = (5.5,0.5) and M = (5.5,5.5) in terms of n/8. (b) The momentum-space 
contour plot of the overlap between \<&i ) and l^f). 



combination of |Oj[) and !<!>£). Without the projection, Eq.rfTOt 
provides a route to construct the BCS ground state and the two- 
particle excited state by using a linear combination of two or- 
thonormal states and l^)- Therefore, we can expect the 
ground state |<I> k ) and the GQP excited state |cD k ) should play 
the same role in the cases under the projection. A further ques- 
tion is whether the GBQP wave function |<t> k Be/) ) mentioned 
above is appropriate to describe the two-particle excitation from 
the (i-RVB state |Oo). In FigQJa), we compare the dispersion 
of four different states discussed above with \<b < ^ BQF ). First, it 
is obvious that |<& k ) exactly reproduces the optimized energy 
of the d-RWB state indicated by the orange dashed line. Based 
on this agreement, we are confident that |<I> k ) should properly 
represent the two-particle GQP excitation. Interestingly, except 
the deviation near the antinodal regions, the energy dispersion 
of |O k ) coincides with \<t>^ B ® P ) very well. We shall return to 
this deviation in the following. Here we see that \®° BQP ) has 
the lower energy than |O k ) around the regions with large d- 
wave pairing amplitude. Therefore, we may conclude that the 
GQP excitation seems reasonable to be constructed by applying 
BQP operators to the d-RVB wave function although |<l> k ) just 
represents the first two-particle GQP excited state. 

Second, we notice that |cD k ) (|<t> k )) shows dispersionless be- 
havior inside (outside) the underlying Fermi surface. This re- 
sults can be easily understood by transforming the representa- 
tion (df) back to (c). In the original representation (c), they are 




Figure 2: The coefficients (a) A k and B^; (b) and D k along the high- 
symmetric momenta. Their definitions are explained in the text. The BCS 
coherence factors and Vj obtained from the optimized parameters are shown 
for comparison. 



given by 



= PNP C clfl H Y] ("q + ^T^J |O><0, 



q^k 



PN r PcY]{ Uq + v q cl/_ ql )\0) {c) . (14) 
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Apparently, there is an additional electron (hole) pair in |O k ) 
(|<£> k )) so that it only can show the dispersion outside (inside) 
the underlying Fermi surface in FigQJa). Next, we study the 
overlap between and |<I> k ) shown in FigQJb). Note that 
without the projection, PFj*) and are orthogonal for every 
momentum. Under the projection the overlap dramatically en- 
hances near the antinodal regions. The shape of the overlap in 
the momentum space is very similar to the li-wave form factor, 
cosk x - cosky, and the maxima about 0.82 right at antinodes. 
It can be easily understood from Eq. (fl~4b in the original rep- 
resentation (c). Since there is an extra Cooper pair c^c^ in 
the overlap, we can expect that Gutzwiller's projection would 
influence the Cooper pair following the of-wave symmetry in 
BCS wave functions. Thus, we can further comprehend the de- 
viation for the excitation energy near the antinodal regions in 
Figllla). 

Furthermore, it is important to examine the details of the 
eigenstates |O k ). We can write down |€> k ) as a linear combina- 



tion of the orthogonal states |<£>. '*) 



(15) 
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where the coefficients and can be easily determined by 
diagonalizing Eq.dTHl. By using Gram-Schmidt method, the 
two orthogonal states |O k '*) can be written as 



<^i<>i^> 



\<&M)\ 2 



(16) 



Thus, the eigenstates |cDj^) are expressed in terms of |0. ^) as 
well, 



(17) 



Here the coefficients Aj^ and are related to and ac- 
cording to Eq.([T6li. Note that the normalization condition guar- 
antees |C^| 2 + \D^\ 2 - 1 but not necessary for Aj^ and B^. 

In order to understand how Gutzwiller's projection affects 
the coherence factors in Eq.([T0l>, we plot the comparison be- 
tween the BCS coherence factors and the coefficients A^ and B^ 
in Figj2fa). We should notice that and \\ are obtained from 
the optimized variational parameters. To avoid the confusion 
arising from the sign of cZ-wave symmetry, we show the abso- 
lute value of vt and A^. Along the nodal direction, interestingly, 
the projection seems not to change the BCS coherence factors. 
This may imply the intimate relation between Gutzwiller's pro- 
jection and the if-wave gap function. However, since the coef- 
ficients AjJ and B^ are not normalized, there is a large enhance- 
ment for the first GQP excited state and small suppression for 
the cZ-RVB ground state around the antinodal parts. In Figj2jb), 
we clearly demonstrate the shape of the coefficients and 
looks similar to the BCS coherence factors in spite of the ex- 
istence of Gutzwiller's projection. The crossing curve of 
and bends more like a hole pocket at {n, jt), and however 
the BCS coherence factor displays a diamond-like underlying 
Fermi surface (not shown). Even so, it is reasonable to believe 
that the GQPs in the <i-RVB wave function are still analogous 
to the BQP picture. 
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4. Conclusion 

In conclusions, we have developed a method to examine the 
idea of BQPs in Gutzwiller-projected BCS wave functions. By 
calculating the two-particle excitation dispersion using VMC 
approach, the GQP excited states have been obtained. We have 
found the GBQP \<&™ QP ) shows almost the same energy as the 
GQP \<&t) except that gives the lower energy around the antin- 
odal regime. The reason is that there exists large overlap be- 
tween and |<I>£) near the antinodes, suggesting the intimate 
connection between Gutzwiller's projection and af-wave Cooper 
pairing. It also results in the deviation of the coefficients in 
|(t)^) from the BCS coherence factors, which might be related to 
the nodal-antinodal dichotomy observed in cuprates by ARPES 
measurements. Therefore, the lower energy around the antin- 
odal regions implies the GBQP excited wave functions are suit- 
able to describe the low-lying excitations in strongly correlated 
superconductors. 
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